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Theoretical background 


haf Introduction 


LJ SupersonicCommercial Transport Aircraft Design 2 cage TI. 
“ Safety ae 4 7 
> Light weight airframe can cause strength, buckling, aeroelastic, and aeroservoelastic issues. a » 
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“* Sonic boom 
> Supersonic flight of “commercial transport” aircraft allowed only over the ocean. 
> Perceived Loudness in decibels 
Y NASA's N+2 goal: 75 PLdB 
Y Concorde: 104 PLdB 
Y High Speed Civil Transport (HSCT): 99 PLdB 
“+ Fuel efficiency 
> Light weight airframe 
> Reduced drag 
L) Developing Quiet Supersonic Technology (QueSST) X-plane 
** Low Boom Flight Demonstrator (LBFD) 
“* Lockheed Martin Skunk Works is the prime contractor for preliminary design. - a 
“+ Loudness: 74 PLdB + QueSST X-plane » 


LJ Major Issue Aa 


** Out-mold-line configuration of an aircraft is design for the desired aerodynamic 
performance. Assume rigid structure. 

** Flexibility of the structure changes the aerodynamic performance. 

** It has been reported that one degree of the tip twist of a LBFD wing and stabilator 
under the cruise flight condition can increase the sonic boom level by 0.2 PLdB and 
1.3 PLAB, respectively. 
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Jig-Shape Optimization Problem Statement 


| Assume unconstrained Optimization 
| Optimization Problem Statement 


?, 


“se Find design variables: {X} = |X, Xo,  Xnay] which 


nsur f 3 


minimize < F(X) = ». AT,* 
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& {AT} = {T}: — {Ta 
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“s {T},= target trim shape at surface GRIDs 


Sonic boom level is computed based on target trim shape. 


“se {T} = trim shape based on design jig shape 
> 


Yigha =e a 
Uigsa = Vigh, + tAjig} 
v {jig}q= design jig-shape 
v {jig}p= baseline jig-shape 
v {Ajig}=jig-shape changes 
{Ajigh = [P]{x} 
v X;= i-th design variable 
Yo [®] = [tbh tode - {Dhnavl 
e {f}; = i-th basis function 


: Eigen vector based on jig shape 
a 


Eigen vectors are normalized as Max deflection = 1 inch. 
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Bs Update Jig-Shape Module: using shape_change.exe 


Lbfd_trim.bdf Sol101.bdf 


From O3 tool variables and basis functions. 


% Yigha = Vigh + [P]{x} 
MSC/NASTRAN “+ Input 
optimization 


Change_trim.exe Sol101.406 
>» design_var: design variables of the current 


Shape.exe optimization step 


Trim.bdf : : ; 
| Grid.dat_— | configuration (a template file) 
Shape_change.exe rine 


Design_var 


i |Basis_functions.dat 


Grid_base.bdf 


Shape.bdf configuration 
Shape.dat 


Grid_update.bdf rd 
Sol103.bdf 


{The 
Differ.exe Target.dat 


{Ta 

(an) 
F(X) 
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= >» AT? 
Intload.dat Nom Oba cele) | 


Extload.dat 


MSC/NASTRAN 


Airload.dat 


Sol103.f06 


J 


j=1 


Structural Dynamics Group Change jig shape using design variables & basis functions. 


L} Shape_change.exe: Change jig shape using design 


> basis_functions.dat: basis functions for the shape 


> grid_base.bdf: GRID information of the baseline 


> grid_update.bdf: GRID information of the updated 
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4) Modal Analysis Module: using MSC/NASTRAN solution 103 


LJ Perform modal analysis using MSC/NASTRAN solution 


From O3 tool Lbfd_trim.bdf Sol101.bdf 103 to change system mass matrix (MGH matrix), 
weight, moment of inertia, and CG location for trim 


LJ Compute six rigid body modes. 


silks bch | Solt01.f06 | 
($}i poli0i-106 


Shape.exe 
{jig}, Grid_base.bdf P 
a, [Griddat__—| 
Shape.bdf 
Shape.dat 


Shape_change.exe 


Grid_update.bdf mi 
Sol103.bdf 


{Te 


Extload.dat Differ.exe 
MSC/NASTRAN (Tha 
Airload.dat ey 
So10506 ner 
= >» AT;* 
| _Sol103.mgh__| Intload.dat To O? tool ie 
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Trim Analysis Module: using ZAERO & change_trim.exe 


L) Change_trim.exe: Update input deck for ZAERO trim 


From O° tool analysis. 
“* Input 
x, Sosmiival MSC/NASTRAN > Lbfd_trim.bdf: template input file file for ZAERO 
= based trim analysis 
Change_trim.exe > Sol103.f06: f06 file from MSC/NASTRAN 
or 3 & Output 


> trim.bdf: updated ZAERO input file to be used for 


7 - Shape.exe 
{jig}p| Grid_base.bdf trim analysis 
Trim.bdf Ts«Grid.dat_— LJ Perform trim analysis using ZAERO 
“+ Input 
Shape_change.exe > Trim.bdf 
Shape.bdf * Output 
Gig}al Grid_update.bdf Seanad > Extload.dat: aerodynamic load + inertial load 


{The 
Differ.exe Target.dat 


Ta 

(ar) 
F(X) 

se 
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= >» AT; 
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Sol103.bdf 


Extload.dat 


MSC/NASTRAN 


i Sol103.f06 _ 


Airload.dat 
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Nok O ka (ove) | 


Intload.dat 
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Objective Function Module: using MSC/NASTRAN solution 101, shape.exe, & differ.exe 


L) Perform static analysis using inertia relief. 


Frat teal (MSC/NASTRAN sol. 101) 
L) Shape.exe: read and write trim results. 
“+ Input 
Xi i MSC/NASTRAN 
pesign var > Trim.f06: MSC/NASTRAN output from sol. 101 
: Change_trim.exe “* Output 
{Phi Grid.dat: GRID geometry information (@ all GRID) 


> 
Shape.exe ‘ Deform.dat: deformed shape (@ all GRID) 


Shape.dat: GRID geom. + deformed shape (@ all 


{jig}, Grid_base.bdf 


rigger GRID) 
> Shape.bdf: shape.dat in MSC/NASTRAN input deck 
Shape_change.exe format 
L) Differ.exe: compute performance index 
., : “+ Input 
ne [—Shape-dat : 
P phabe:dat > Shape.dat 


{T}> Target.dat: {T}, @ surface GRID 


Differ.exe Target.dat 


Ma * Output 
> Pindex.dat: performance index for objective 
GO. 


function; F(X) 
by : Tvect.dat: {AT} 
Intload.dat To O? tool Z 


Sol103.bdf 


MSC/NASTRAN 


Extload.dat 


Airload.dat 


Sol103.f06 


> 
AT)", Shape_diff.dis: {AT} @ surface GRID for 
MSC/PATRAN plotting 
> Surface.dat: shape.dat @ surface GRID; {T}q 


~ Compute trim deformation using MSC/NASTRAN solution 101 using inertia relief. _ Chan-gi Pak-9 
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Gy Compute Starting Design Variables: Using Least Squares Surface Fitting Technique 


O {AT}, = {7}. — {Th 
“* {T},= target trim shape at surface GRIDs 


“* {T},= trim shape based on the baseline jig-shape 


trim analysis 
UI Fitting {AT}, surface using perturbed shapes {AT};, i = 1, 2,...,ndv 


?, 


“* Perturb baseline jig-shape using basis functions [®] 
> Yigha = Vighy + [P]{X} 
> Where, {¢},= i-th basis function 
PUD TAO) 
> {AT}, = {T}, — {T}, (i-th perturbed shape) 
“* Define a matrix: [W] = [{AT}, {AT}, ... {AT} nav 
O [P]{X} = {AT}, 
oe [W)"[P]{X} = [VY] {AT}, 
(Pe) Te Ped = (TP) Py 473, 
QO) Starting design variables: {X} = ([W]" [W])~7[W]? {AT}, 
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Computational validation 
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Structural Finite Element Model 
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Summary of Natural Frequencies (Baseline Configuration) 


Frequency | Frequency (Hz) 


Notes 
| 5.634 634 | | First fuselage bending 
= 9.045 pp First wing bending + forward fuselage vertical bending + stabilator rotation 
11.97 Lea Forward fuselage vertical bending + first wing bending + stabilator rotation (Asymmetric) 
pa 19.23 a a Wing tip bending + T-tail rotation + flap bending (Asymmetric) 


ce 20.08 | | ttaitrotation(Asymmetric) T-tail rotation (Asymmetric) 


20.54 Wing tip bending + T-tail rotation + aileron rotation + flap bending + forward fuselage vertical bending 
Asymmetric 
21,75 a a Aileron rotation + flap rotation + T-tail bending + outboard wing bending torsion 
= 22:16 a Flap rotation + aileron rotation + wing tip bending + T-tail bending (Asymmetric) 
22.70 a Flap rotation + aileron rotation + T-tail bending (Asymmetric) 
296 | | | Tealtaending (Aaymmetri 
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Mode 7: 5.634 Hz Mode 9: 9.045 Hz 


1.294 1.468 
1.207 1.370 
1.121 1.272 
1.035 1.174 
0.949 1.076 
0.863 0.978 
0.776 0.881 
0.690 0.783 
0.604 0.685 
0.518 0.587) 
0.432 0.489 
0.345 0.392 
0.259 0.294 
0.173 0.196 
0.087 0.098 
4.90-004 2.74-004 


Symmetric first wing bending + forward fuselage vertical bending + 
horizontal tail rotation (in-phase: forward fuselage & wing) (out-phase: 
wing and horizontal tail) Chan-gi Pak-14 


first fuselage vertical bending 


Mode 11: 11.97 Hz Mode 15: 14.76 Hz 


1.389 3.248 
Asymm etric 1.296 3.031 
1.204 2.815 
1.111 2.598 
1.019 2.382 
0.926 2.165 
0.834 1.949 
0.741 1.732 
0.648 1.516 
0.556 1.299 
0.463 1.083 
0.371 0.866 
0.278 0.650 
0.186 0.433 
0.093 0.217 
7.79-004 2.49-004 


Symmetric first wing bending + forward fuselage vertical bending + 
horizontal tail rotation (out-phase: forward fuselage & wing) (in-phase: 
wing and horizontal tail) 


Symmetric horizontal tail rotation 
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Mode 17: 19.23 Hz 


Asymmetric 


symmetric wing tip bending+Ttail rotation + flap 


Asymmetric nian 


Mode 19: 20.08 Hz 


2.688 
2.509 
2.330 
2.151 
1.972 
1.793 
1.614 
1.435 
1.256 
1.076 
0.897 
0.718 
0.539 
0.360 


0.181 


symmetric ttail rotation 


6.808 


6.354 


5.900 


5.446 


4.992 


4.539 


4.085 


3.631 


3.177) 


2.423 


2.270 


1.816 


1.362 


0.908 


0.454 
3.67-004 
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“3 Mode 20: 20.54 Hz Mode 22: 21.75 Hz 


Asymmetric 5199 


2.055 
1.918 2.052 
1.781 1.906 
1.644 1.759 
1.507 1.613 
1.370 1.466 
1.233 1.320 
1.096 1.173 
0.959 1.026 
Asymmetric 
0.822 0.880 
0.685 0.733 
0.548 0.587 
0.411 0.440 
0.274 0.294 
0.137 0.147 
3.86-004 6.63-004 


symmetric wing tip bending + ttail rotation flap & airleron rotation + 
forward fuselage bending + nose landing gear vertical bending (out-phase symmetric airleron + flaperon (in-phase)+ttail(pitch +yaw) 
wing tip & forward fuselage) (out phase wing tip & ttail) . Sisco 


“3 Mode 23: 22.16 Hz Mode 25: 22.70 Hz 


Asymmetric ane 3.299 
2.299 3.079 
2.135 2.859 
1.971 2.639 
1.807 2.419 
1.642 2.199 
Asymmetric 

1.478 4 1.979 
1.314 1.759 
1.150 1.540 
0.986 1.320 
0.822 1.100 
0.658 0.880 
0.493 0.660 
0.329 0.440 
0.165 0.220 
9.25-004 4.92-004 

symmetric Flaperon+airleron (out-phase) +ttail(pitch+yaw) +forward symmetric flaperon+airleron (out-phase)+ttail(pitch+yaw) + forward 

fulage and airleron(in-phase) fulage and airleron(out-phase) 
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Nf Trim Shape Difference (Baseline Configuration) 
LAT to 


L} Weight: 
“* Cruise = 18499.99 lbf 
WU) Forward CG location 
** x=836.09 inch, y=-0.1897 inch, z=100.68 inch 
Mach: 1.42 
Altitude: 55000 ft 
Aileron deflection angle: 0.5 deg 
T-tail deflection angle: 6.47 deg 
{AT}.o = {7}. —{T}o 
“* {T},= target trim shape at surface GRIDs 


DOUOCD 


“* {T},= trim shape based on optimum jig shape 
Y {figho = Vigh + [P]{Xto 


trim analysis 
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Asymmetric 
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Use least-squares surface fittin 


OBJ = 17.48 
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Start Configuration 
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inch 
0.240 


0.224 


0.208 


0.192 


0.176 


0.160 


0.144 


0.128 


0.112 


0.096 


0.080 


0.064 


0.048 


0.032 


0.016 
6.42-005 


GY Optimization #1:{AT},, = {7}, —{T}o 


Use Optimization 


OBJ = 17.34 


O 


timum Configuration 


Mode 37: 30.79 Hz Mode 48: 42.96 Hz 


6.864 Asymmetric 8.299 
6.406 7.745 
5.949 7.192 
5.491 6.639 
5.034 6.086 
4.576 5.532 
4.118 4.979 
3.661 4.426 
3.203 3.873 
2.746 3.320 
2.288 2.766 
1.831 2.213 
1.373 1.660 
0.916 1.107 
0.458 0.553 
3.90-004 1.92-004 
symmetric canard bending symmetric ttail 
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Optimization #2: {AT},,= {T}. — {T}, 


Use least-squares surface fittin 


OBJ=6.148 inch 


0002 
0.086 
0.080 
0.074 
| 0.068 
rs [osare - 
0.055 
ra [nase 

0.049 
0.043 


0.037 11 


12 


os 


0.031 


ea 


0.025 
0.019 
0.012 


0.006 
Start Configuration 1.43-004 


Use Optimization 
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inch 
0.085 


0.079 


0.074 


0.068 


0.062 


0.057 


0.051 


0.045 


0.040 


0.034 


0.028 


0.023 


0.017 


0.011 


0.006 
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GY optimization #3: {AT},.= {T}; — {T}o 


Use least-squares surface fittin inch 
OBJ=0.09566 
0.0064 
0.0059 
0.0055 
0.0051 
0.0047 
| 6 | 0.6221 | 0.0043 
0.0038 
8 | -.0265 | 0.0034 
5 [000167 
0.0030 
0.0026 
0.0021 
0.085 inch 0.0017 
0.0013 
0.0009 
0.0004 


Start Configuration 
4543-006 


Use Optimization 


.O 
.O 
0.1165 
-0 
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Optimization Results 
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Nf Optimum Aircraft Configuration 
{AT}, inch 


“L20.152” Weight (Ib) | 18499.99 | 18499.97 | 0.00 | 
NESE X-C.G. (inch) | 836.0883 | 836.0881 | 0.00 ‘| 


[in [629919800.] 629901800,[ 000 


R: -0.154" SO. 
L: -0.151” \ 


: Baseline Configuration > 
: Optimum Configuration I] 
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Summary of Natural Frequencies before and after optimization 


Frequency (Hz) 
% Notes 
5.634 5.633 First fuselage bending 
po 9.045 9.034 First wing bending + forward fuselage vertical bending + stabilator rotation 
11.97 11.97 | 0.00 | Forward fuselage vertical bending + first wing bending + stabilator rotation (Asymmetric) 
19.23 19.23 | 0.00 | Wing tip bending + T-tail rotation + flap bending (Asymmetric) 


20.08 20.08 | 0.00 | T-tail rotation (Asymmetric) 


20.54 20.55 Wing tip bending + T-tail rotation + aileron rotation + flap bending + forward fuselage vertical bending 
Asymmetric 

21,75 21.76 Aileron rotation + flap rotation + T-tail bending + outboard wing bending torsion 

22:16 2247 Flap rotation + aileron rotation + wing tip bending + T-tail bending (Asymmetric) 

22.70 22.70 | 0.00 | Flap rotation + aileron rotation + T-tail bending (Asymmetric) 

30.79 30.76 Canard bending 

42.96 42.97 T-tail bending (Asymmetric) 
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Conclusion 


UI In this study, the jig-shape optimization is performed using the two step approach. 
“* The first step is computing the starting design variables using the least squares surface fitting technique. 
“* The next step is the fine tune of the jig-shape using the numerical optimization procedure. 
** Assume unconstrained optimization 
> The maximum frequency change due to the jig-shape optimization is less than 0.12%. 


> The minor changes in mass moment of inertia are observed. (mostly less than 0.07%; maximum 0.73%) 


LJ Thirteen basis function are used in this jig-shape optimization study. 
“* Total of twelve symmetric mode shapes of the cruise weight configuration. (Asymmetric shapes exist) 


** A residual shape is also selected as a basis function. 


LJ The maximum trim shape error of 1.174” at the starting configuration becomes 0.0088” at the end of the third optimization run. 
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